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ABSTRACT 

Equations of motion have been derived for a fin-stabilized 
projectile with time as independent variable. They have been 
simplified assuming that the rates of yaw, roll and pitch, and 
also angle of attack, sideslip angle and angle of roll to be 
small. The aerodynamic forces and moments are represented by 
the corresponding coefficients. These coefficients have been 
estimated for a 120 mm i artillery projectile and the equations 
have been solved using a variable-order, variable-step Gear 
method. Corrections to the trajectory due to effects of atmos- 
pheric variations and wind have been incorporated. Parametric 
study consisting of variations in aerodynamic coefficients has 
been carried out to study their effect on the trajectory. The 
results have been compared with those in the range table. 



CHAPTER 1 

INTRODUCTION 


1.1 GENERAL 

Rapid advancement in technology since World War II has 
made the present day battle field environment highly complex. 
The modern armies are equipped with weapon systems which are 
both accurate and lethal and their increased mobility can 
greatly enhance the dimensions of the battle zone. Future 
conflicts will be short, fast and fluid. The lethal weapons 
which are the potential targets of artillery will be deployed 
swiftly and moved to a new location before they can be fired 
upon. Fleeting opportunities offered by the enemy will have 
to be seized to destroy the targets. Speed will be vital in 

this regard. Neutralization of targets by artillery will 

* 

impose exacting demands on the accuracy of target acquisition 
devices, engagement procedures and weapon systems. This 
necessitates the employment of an efficient command, control, 
and communication system which can provide sufficient target 
intelligence, rapid target processing including assessment of 
priorities, optimum allocation of artillery resources, accurate 
computation of gun data, and reliable and secure communica- 
tion, all with high speed of response. Only then it will be 
possible to use artillery in its vital battle winning and 
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challenging role. 

The factors which greatly influence the accuracy of 
artillery fire are : 

. The accuracy of the location of gun and target 

. The accuracy, reliability and speed of communication of 

the firing data 

. The accuracy of the meteorological report which includes 
information about the pressure, temperature, density and 
wind speed and direction in various layers of the 
atmosphere 

. The accuracy of the weapon system 
. The training standard of the crew. 

1 .2 MOTIVATION 

With the help of survey techniques and target acquisition 
devices like radars, sound ranging, laser range finders and 
remotely piloted vehicles, the gun and target locations can be 
determined with high degree of precision. Meteorological radars 
equipped with radio-sonde and on-line computer to process the 
radio-sonde data will soon replace the balloon theodolites to 
provide accurate and timely information about the weather 
parameters. To achieve first salvo-effectiveness of artillery 
fire, it therefore becomes imperative to calculate the firing 
data for the guns quickly and accurately. 
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The present method of determining the firing data is slow 
and prcne to human errors. Range table is used to compute the 
angle of projection for the guns to achieve the desired range. 
There is a range table supplied with each type of the artillery 
gun by the manufacturer. These range tables are compiled by 
actually observing the fall of shot at certain standard eleva- 
tions. For intermediate elevations the data are interpolated. 
Corrections for the non-rigidity of the trajectory and other non- 
standard conditions are also included in the range table. 

In a fluid battle situation where target information flows 
in continuously, paging through the range table to find the 
firing data for large number of targets may considerably slow 
down the response of artillery guns. A computer can greatly 
improve the accuracy and speed of response in this field. Range 
table contains large volume of data and since there is a range 
table for each type of projectile it is uneconomical to store 
the data of all the range tables in a computer memory and use 
it to calculate the firing data. Using a computer merely as a 
substitute for range tables will definitely improve the speed 
of response. But the disadvantages inherent in compilation of 
range tables and large requirement of memory space prohibit 
this approach. The most effective use of a computer lies in its 
ability to compute the firing data for all types of guns for 
different elevations and muzzle velocities. This can be achieved 
only when it is programmed to solve the equations of motion and 
determine the actual trajectory of the projectile through the 
atmosphere. 
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1.3 CURRENT APPROACH 

In the present investigation a technique analogous to 
that used for the study of aircraft flight has been employed 

to derive the equations 0 f motion of a pro jectile fired from 

artillery gun. The anal ysis is basically describing the motion 

of the projectile in terms ^ , 

ms of a set of differential equations 

written in an inertial f v. am „ ^ - . . . ' 

r rame of reference under the action of 

external forces and moment,. t +. u , . 

nt s namely the aerodynamic and gravity. 

The aerodynamic forces and moments acting on the projectile 
are represented in terms of the corresponding aerodynamic 

coefficients. These equati ^v,,. +.i_ . , .. 

H^arions are then suitably linearized 

and non-dimensionalized. 

This method of computing the trajectory from the equations 
of motion is a more genera! approach and can be applied to any 
artillery projectile with known aerodynamic characteristics. It 
offers the following advantages : 

The firing data can be computed quickly and reliably 
. Complete information about the trajectory of the projectile 
can be known. Any modifications like inclusion of thrust 
effect etc. can be easily taken care off. 

It is more economical than the present method of predicting 
the trajectory which entails maintaining establishments to 
compile the range table by actually observing the fall of 
shot. The process has to be repeated for every new 
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projectile which differs in shape and size. However, 
with the new method this expensive field study can be 
eliminated. Only aerodynamic coefficients need to be 
determined in a laboratory or estimated analytically. 
Parametric study can be carried out to arrive at the 
optimum combination of muzzle velocity and angle of 
projection for engaging a target. 

The parametric study may also help in arriving at the 
optimal projectile shape and inertial properties. 
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CHAPTER 2 

BASIC CONCEPTS 

In this chapter the geometry of fin- stabilized projectile 
and its orientation in space are introduced and the conditions 
of flight are postulated. The two axis systems, the earth- 
fixed axis system and the stability axis system, required for 
the development of equations of motion are discussed. The 
transformation matrix relating the two axis systems and the 
equations relating the components of angular velocity of 
projectile to the Eulerian angles of projectile in space have 
also been included. 

2.1 GEOMETRY OF THE PROJECTILE 

The projectiles fired from artillery mortars have 
stream-lined shape provided with fins at the tail end to 
provide stability during flight. These projectiles are 
considered to have a 90° roll symmetry. The geometry of a 
typical artillery projectile and its orientation in space are 
shown in Fig. 2.1. Consequence of this 90° roll symmetry is 
to reduce the number of aerodynamic coefficients in the equa- 
tions of motion and thus reducing their complexity. 

2.2 COORDINATE SYSTEM 

In developing the equations of motion two right-handed 
orthogonal coordinate systems have been chosen. One of them is 
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stability axis system which is fixed to the projectile. The 
other set is fixed with respect to the earth and is known as 
earth-fixed axis system. If the rotation of earth is neglected 
the earth-fixed axis system forms an inertial system of 
reference frame. To derive the equations of motion Newton's 
laws are applied in this space fixed coordinate frame. The 
physical quantities such as linear and angular velocities and 
aerodynamic forces etc. are conveniently dealt with in a 
coordinate frame moving with the body and fixed to it. The 
equations of motion written with reference to space fixed 
coordinate system are subsequently transformed to stability 
axis system. Knowing the relative orientations of the two 
axis systems, the variables of interest in one coordinate 
system can be transformed to another coordinate system. The 
Fig. 2.1 shows the two axis systems and their relative orien- 
tations. 

2.2.1 Space Fixed Axis System 

It is designated by X 0 Y 0 Z q and has its origin at the 
point of launch of the projectile. The directions of the 
individual axes are as follows : 

X Q : aligned parallel to the longitudinal axis of the 

projectile and is positive in the forward direction. 

Y 0 : is horizontal at 90° to X 0 and is positive to the 

right of an observer facing forward. 

Z Q : points downwards perpendicular to X Q and Y Q . 
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2.2.2 Stability Axis System 

It is fixed to the projectile and is designated by XYZ« 

The origin of this axis system coincides with the centre of 
gravity of the projectile and the directions of the individual 
axes are as follows : 

X : coincides with the longitudinal axis of the projectile 
and is positive in the forward direction. 

Y : is at right angles to X and initially points in the 
same direction as Y Q . 

Z : is at right angles to X and Y and initially, points in 

the same direction as Z . 

o 

By virtue of the 90° rotational symmetry of the projectile 
all the three axes of this axis system are the principal axes 
of inertia. 

2.3 ORIENTATION OF THE PROJECTILE IN SPACE 

The orientation of the projectile in space at any instant 
of time is defined by three angular coordinates ¥ , 9 and 0 as 

t 

shown in Fig. 2.1. These angles called Eulerian angles des- 
cribe the orientation of stability axes, XYZ f with respect to 
space-fixed axes, X 0 Y 0 Z Q . They are zero initially when the 
two coordinate systems coincide with each other and any sub- 
sequent orientation of the projectile can be obtained by 
giving three consecutive rotations of magnitudes S , © and 
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in order. Thefce rotations are non-commutative and must be 
carried out in an order always. The order generally followed 
is as given below : 

. Starting from XYZ coincident with X o Y Q Z 0 rotate XYZ 
through an angle ¥ about the Z Q axis. This results in 

x lYl z r 

• Rotate X^Y^Z^ about Y^ axis through an angle ©. This gives 

* 2 * 2^2 ^2 position. 

. The third rotation is about X 2 (X) axis by an angle 0. 

This brings the Y and Z axes to their final positions and 
the projectile to its actual orientation in space. 

The Eulerian angles if , 9 and 0 can also be visualized as 
follows : 

If : angle between X Q axis and the line of intersection of the 
planes X Q Y 0 and XZ Q . 

9 : angle between the line of intersection of the planes X Q Y 0 

and XZ Q and X axis. 

0 : angle between the line of intersection of the planes YZ 

and X q Y 0 and Y axis. 

Figure 2.1 shows the direction of gravity vector which is 
vertical. Since Y Q axis is horizontal the X Q Z Q plane becomes 
the vertical plane and hence contains the gravity vector. The 
X Q axis is inclined to the horizontal plane at an angle © 0 ,the 
angle of projection of the projectile. The angle between Z Q 
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axis and the gravity vector is also equal to angle ©o. This 
angle is considered positive when X 0 axis is above the hori- 
zontal plane. In case of artillery projectiles angle 0 Q will 
always be positive. 

2.4 LINEAR AND ANGULAR VELOCITIES 

The orientation of the velocity vector V with respect 
to stability axes X,Y,Z is given by the angle of attack a and 
the angle of side slip p as shown in Fig. 2.2. The angles a 
and p are defined by 


a = 90° - * V Z, 
p - 90° - * V Y . 

From these definitions of angles a and p it follows that 
the components of velocity V along X,Y and Z axes are 


V X - v 


J l-sin^a-sin^p , 


and 


/ Y = V sin p, 


Vj, = V sin a . 


( 2 . 1 ) 


The angular velocity components p»q»r of the projectile along 
X,Y f Z axes are related to the Eulerian angles f © and 0 , and 
their time derivatives ¥ , 0 and 0 by the following kinematic 
relations [1] ! 



u 


p * $-■- W sin© 

• , • . , ( 2 . 2 ) 

q = © cos0 + 't cos© sm0 

r = * cos© cosjZf - 6sin0 

2.5 CONDITIONS OF FLIGHT 

Following conditions are postulated for the geometry of 
the projectile and its flight : 

. The projectile has 90° rotational symmetry about X axis 

. The angles and 0 and also the time derivatives 

* * • • 

a, {3, Sf , © and 0 are small. 

It follows from the second assumption and from eqn. (2.2) that 
the components of angular velocity, p,q and r are also small. 

2.6 RELATIONSHIP BETWEEN THE TWO AXIS SYSTEMS 

The equations of motion are written in XYZ axis system. 

But in order to analyse the trajectory the distances travelled 
by the projectile must refer to X 0 Y Q Z 0 axis system. The trans 
formation from one axis system to another can be made using the 
relationship between the unit vectors of the two axis systems. 
This relationship between i,j,k and i 0 >j 0 »k 0 is given by [l] 
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— * 

*0 


cos© cosy 

cosy sin0 sin© 
-siny cos0 

cosy cos0 sin© 

+s iny sin0 


i 



cos© sin 4 

siny sin0 sin© 
+cosy COS0 

siny cos0 sin© 
-cosy sin0 


j 



-sin© 

sin0 cos© 

cos© COS0 


k 






(2.3) 


where i,j,k and i o ,j o ,k Q are unit vectors in XYZ and ^ 0 Y Q Z 0 axis 
systems respectively. If angle 0 is small as postulated in 
Sec. 2.5, sin0 can be approximated by 0 and cos0 by 1. With 
this eqn. (2.3) reduces to 


L 

i o- 


cos© cosy 0cosy sin© - siny cosy sin©+0siny 


i~ 


SB 

cos© siny 0siny sin© + cosy siny sin©-0cosy 


j 

_ k o_ 


-sin© 0cos© cos© 


k 


(2.4) 
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CHAPTER 3 
EQUATIONS OF MOTION 


In this chapter the general equations of motion of a 
projectile are derived with time as independent variable. 
After suitable transformation, linearization and nondimen- 
sionalization these equations are rearranged in a form which 
is convenient for trajectory analysis. 

3.1 THE VECTOR EQUATIONS 


Newton's second law of motion expressing the conservation 
of linear and angular momenta provides the two basic vector 
equations which represent the motion of the projectile. They 
are 


dt r R • 

and 

Ht » = % * 


(3.1) 

( 3*2) 


where f N "iand H represent linear and angular momenta respectively, 

mmm 1 

the resultant external force acting on the projectile includ- 
ing the thurst and the resultant external moment. 

The components of linear momentum f^along the X,Y and Z 
axes are mV^, mVy and mV^ respectively so that 


imV^ + jmVy + kmV z , 


(3.3) 
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where m is the mass of the projectile and V^, Vy and V z are 
the velocity components given by eqn. (2.1). 

Similarly the angular momentum vector H can be written 

and hL, as 


in terms of its components Hy 


H * iH x + jHy + kH z , 


(3.4) 


where 


p I X “ rl xz “ qI XY» 

Hy — qly “ Plj^y “ rXy-ry 


(3.5) 


^ ”* rI z - ql yz - pl xz ♦ 

In eqn. (3.5) I^IyjI^ are the moments of inertia of the 
projectile about X,Y,Z axes and I^yt I xz and Iy Z are the 
products of inertia. By virtue of the 90° rotational symmetry 
Iy and I z will be equal and all the products of inertia will 
be zero. Thus the components of H given in eqn. (3.5) reduce 


H y - ply» 


^y * qly » 


(3.6) 


Hz - rl z . 

And hence eqn. (3.4) reduces to 

H » ilyP + jl v q + kl 7 r, 


(3.7) 
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where p,q and r, the components of angular velocity are given 
by eqn. (2.2) • 

The unit vectors in eqns. (3.3) and (3.7) are functions 
of time and hence the time derivatives in eqns. (3.1) and (3.2) 
have contributions arising from this which are accounted for by 
using the Poisson’ s formula 


dA — — — 
^pjr — A + a) x A 


(3.8) 


where, 


i any vector quantity with A^Ay and as 

components along X,Y and Z axes respectively, 

: rate of change of A with respect to X Y Z, the 
rotating axis system. 


and 


dA 




rate of change of A with respect to X 0 Y o Z Q , 
the fixed axis system, 

angular velocity of the rotating axes X Y Z. 


In the component form eqn. (3.8) can be expressed as 


d Ay 

a** 


« A 


x + qA z - rAy , 


dAv 

eft 


Ay + rAj£ - pA z , 


(3.9) 


dA 


3T Z = + pAy - qA x 


and 
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3.2 FORCES AND MOMENTS 

The external force F R in eqn. (3 f !) comprises in general 
the aerodynamic force F^, gravitational force and the thrust 
force Fj. so that 

Fr + T G + h < 3 - 10 > 

The projectile under consideration has no thrust force acting 
on it and, therefore, Ej. in eqn. (3.10) will be taken to be 
zero. The resultant external moment in eqn. (3.2) is only 

aerodynamic moment. 


3.2.1 Aerodynamic Forces and Moments 

The aerodynamic forces and moments can be expressed in a 
general form as 


A X,Y,Z 


9x,y>z 


v*s 


(3.11) 


whore suffixes X,Y,Z refer to the components of the force F^ or 
the coefficient C in the X,Y,Z directions. Similarly the aero- 
dynamic moments are expressed by 

m X,Y,Z = C l,m,n 4 ^ (3 - 12) 

where l,m,n refer to the components about X,Y,Z axes respecti- 
vely. 
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In eqns. (3.11) and (3.12) is the density of air and 

S and L denote the characteristic area and the characteristic 

length of the projectile respectively. Coefficients Y,Z anc * 

C, are in general, functions of the variables a,p ,a,f3 ,p,q 

x f m,n * 

and r and can be expanded in Taylor's series in these varia- 
bles. A general expression for C Q can thus be written as 


C = C 
a a, 


+ C a+ C 

cln 

a p 


6 + c 
K ^ ♦ 


2V + u a£ 2V * 


- Ek + r Elk 

a p 2V \ 2V 


+ C 5k [higher order terms] (3.13) 

a r 

where ’a’ stands for X,Y,Z,l,m,n . A quasi-linearity approxi 
mation is made so that the coefficients in the series are 
independent of the variables. These coefficients are called 
aerodynamic derivatives or stability derivatives and are in 
general, functions of the Reynolds number and Mach number. The 
stability derivatives are defined as follows : 
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It can be seen from eqn. (3.13) that each force or moment 
coefficient requires 8 stability derivatives for representa- 
tion. In order to describe all the force and moment coeffi- 
cients completely, 48 aerodynamic derivatives will be required. 
By virtue of 90° rotational symmetry of the projectile this 
number will reduce considerably and this simplification is . 
discussed in the next paragraph. 

3*2.2 Consequences, of Mirror and 90° Rotational Symmetry 


Consider two motions of the projectile so that one is the 
mirror image of the other about the plane of mirror symmetry-XZ* 
Then the aerodynamic force and moment due to the second motion 
must be the mirror images of those due to the first motion. To * 
illustrate, consider the term rolling moment due to 

rolling velocity, and C , the pitching moment due to 

m p zv 

rolling velocity. The mirror image of rolling velocity p is 

always - p. The mirror image of C. (§k) should be - C, (§77) . 

P dL? 

This is indeed so since changing p to -p changes C^ (5^) to 

P 

- Cj (57)* The stability derivative C^ is therefore allowable 
and can have a non-zero value. On the other hand the mirror 


image of the pitching moment M y is M y , hence the mirror image of 

C ra <%) should be C m (|t) . it can be seen that this term 
P P 

actually reverses its sign when p is replaced by its mirror 
image -p.This clearly violates the requirement that due to 
mirror symmetry the mirror image of pitching moment My should 
be My itself* Therefore the aerodynamic derivative C m is not. 
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allowable and must be equtql to zero,. By applying the criterion 
of mirror symmetry similarly to all the stability derivatives 
the following derivatives are found to be equal to zero identi- 
cally : 



Cy = Gy = Cy 9 = Q, — 0 
o *a l a q 





(3.14) 


90° rotational symmetry implies that a rotation through 90° 
transforms the projectile into itself. It is obvious that the 
aerodynamic derivatives must have values such that a rotation of 
vector V and u> by 90° about X axis will produce similar rota- 
tion of vector F R and-M R by 90° about the same axis. 


Let A be any vector and A* be the vector obtained by 
rotating A by 90° about the X axis (ref. Fig. 3.1). If the 
direction of rotation is taken from Y to Z as shown in the 
figure, then the X component A x does not rotate, the component 



20 


Ay moves to Z direction and the component A z moves to -Y 
direction so that the new components Ay , Ay , A 7 are related 

I 4 ^. 

A^r > Ay * 


A x* ~ A x * 


a y*. a z * 


(3.15) 


and 




The relationship between the components of the aerodynamic force 
and moment and also those of V and w after rotation of V and w 
can thus be obtained by applying eqn. (3.15) and they are 

(V*, w^.) = F^( V,w) f 
Fy ( V* » ) — **F z (V»to) f 


Fy <V*,W*) = Fy(V,W) , 


(3.16a) 


(Vj(.»m^.) — M^( V,u)) f 
My (V*,W*) = -M z (V,w) , 

V 

(V*, <*>#.) = My(V # (d) , 
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and 


V v 


V, 


cl 


- P 




= -a 


P* = P 
q» = -r 


(3.16b) 


It is noticed that the values of the force and moment components 
after rotation must be computed using the rotated values of 
velocities V* and w* as indicated by the functional notation of 
eqn. (3.16a). Equation (3.16a) must hold for all values of 

V and w. Hence by computing the components of forces and 
moments from eqn. (3.16a), substituting for the components of 

V and to from eqn. (3.16b) and then using eqns. (3.11) - (3.13) 
the following relationships for the stability derivatives are 
obtained : 




* 2 . 


a 




~ C Z* 

a 




=, a 


c = -c 


n 


m_ 


n f3 


-C m . 

m a 


m 


(3.17a) 
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C X “ ~ C x “ ~ C Z ~ C l p ~ C 1 ~ C 1 

a a q p o p p r 


: = C =0 

m n 
o p 


(3.17b) 


Thus the plane of symmetry assumption results in reducing all 
the coupling derivatives between the longitudinal and lateral 
modes to zero as given in eqn. (3.14). The 90° rotational 
symmetry further simplifies the problem by establishing equa- 
lity between some of the longitudinal derivatives and lateral 
derivatives and also by reducing some of them to zero as given 
in eqn. (3.17). 


Finally, the components of aerodynamic force and moment 
can be expressed as follows : 


F A x = 2 C X c 


% = I v2s(< Vp P ' + 27 + C Y r f? 5 

F a = 2 a + C 7 , + C 7 ^) 


J Z' 2V " 2V 
a q 




V^SL C x 


Ek 

2V 


(3.18) 


My 


3 V 231 - < c m « + c m- § + C m # 


M z = | ^SUC^fS + C„. § + 


In view of eqn. (3.17a) there are only 8 independent aerodyna- 
mic derivatives in eqn. (3.18). 
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3*2.3 Gravitational Force 

The components of gravitational force along X Q , Y Q and Z Q 
axes can readily be determined from Fig. 2.1. It can be seen 
that Fq lies in the X Q Z 0 plane and that Y Q is horizontal. 
Therefore, 

F q » i o (-mg sin0 Q ) + k Q (mg cos0 Q ) . (3.19) 

The components of Fq along X,Y and Z axes can now be found by 
using the rotation matrix of eqn . (2.4) and are given by 

Fq = -mg(sin9 Q cos9 cos $ + cos© Q sin0) , 


Fq^ = mg[0cos© Q cos9 - sin© 0 sin9 cos ¥ ) + sin0 Q sin^f ], 


and 


F n = mg[cos© cos© - sin© (cos ¥ sin© + 0sinS)] 


(3.20) 


3.3 THE DYNAMICAL EQUATIONS 


The equations of motion are derived from the basic vector 
eqns. (3.1) and (3.2) by resolving the vector along the coordi- 
nate axes and using the equations for the external forces and 
moments. The final six scalar equations are obtained following 
the steps outlined below : 

. Equations (3.1) and (3.2) are transformed using eqn. (3.9) 

. Components of linear and angular momenta are introduced 
from eqns. (3.3) and (3.7) 
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• The linear and angular velocity components are substituted 
from eqns. (2.1) and (2.2) in terms of V,a,p,¥ ,9,0 and 

‘ », 0 , 0 . 

. The components of aerodynamic force and moment from eqn. 
(3.18) and those of gravitational force from eqn. (3.20) 
are finally substituted. 

3.4 LINEARIZATION 

It is assumed that the angles, 0, a, p and the derivatives 

• * * ♦ 

a, , 9 , and 0 are small and hence only the first order terms 

in these variables are retained and the higher order terms are 
neglected. The procedure followed, therefore, is to group the 
first order and higher order terms separately in each equation 
and then neglect the higher order terms in comparison to the 
first order terms. 

3.5 NON DIMENSION ALIZATION 

The final step in the derivation is to non-dimensionalize 
the equations using the muzzle velocity and the charactoristl* . 
dimensions of the projectile. The non-dimensional time t + is 

i 

defined as t = 7 , where t is the characteristic time given 

V 

by x = , V Q being the initial velocity of projection of the 

projectile and L its characteristic length, t indicates the 

time of flight for one characteristic length of the projectile. 

The non (dimensional speed U is defined as U = and the non- 

d° 

dimensional time derivative operator D as D = — . 
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The details of the derivation, linearization and non- 
dimens ionalizati on of the equations are given in Appendix A and 
their final forms are as follows : 

. X Force Equation 

DU . - U 2 + F(sin© 0 cos© cos SI + cos© 0 sin©) = 0 (3.21) 

• Y Force Equation 

(i - ^ c,.) UDP + cosed - $jj c^Judi + 

- F(sin© 0 cos© cos^f + cos© Q sin©)p - 

F[0(cos© o cos© sin© Q sin© cos¥ ) + sin© Q sinS ] = 0 

(3.22) 

• Z Force Equation 

k c z' a )uoa - ll+ k %> UDS + k (c >c 0 - c z a )u2 “ 

- F(sin© 0 cos© cos W + cos© Q sin©) a - 

F.[cos© 0 cos© - sin© 0 (cosf sin© + 0 sinSf )] = 0 (3.23) 

. X Moment Equation 

D 2 0 - sin© D 2 f - K,C, UD0 + K, C, sin© UD* = 0 (3.24) 

P P 

• Y Moment Equation 

D 2 © - IUC-. UDoc - KgC UD© - 2^^ l^a = 0 (3.25) 

a q a 
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K^Ch cose UD¥ - 21^0 U 2 ^ = 0 

P (3.26) 


„ - PSL° _ PSL 

^2 “ 4I y “ 41 z * 

3.6 THE KINEMATIC RELATIONS 

Equations (3.21) - (3.26) are six differential equations 
in six unknowns U,a,p,¥ ,9 and 0 and a solution of these gives 
the velocity components along X,Y and Z axes. However, for the 
description of the trajectory of the projectile we need to 
transform these components along X 0 Y q Z 0 system so that they can 
be integrated with time in the space-fixed X Q Y 0 Z 0 frame of 
reference to get the trajectory. Using transformation matrix of 
eqn • (2.4) the relationship between the velocity components in 
XYZ and X Q Y 0 Z 0 axis systems can be obtained as follows : 


Z Moment Equation 

cos© D 2 ¥ -L C . UDp 

n p 


where 


F = , 

V o 


m_ 

PSL » 


K, 


PSL J 

41 * 
X 
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”*o 


cos© cos * 

0cos * sin© - sin* 

cos* 

sin© + 0s in * 


"Vx 

Xo 


cos© sin* 

0sin* sin© + cos* 

sin* 

sin© - 0cosW 


V Y 

A 


-sin© 

0cos© 

cos© 



_ V Z_ 


(3.27) 

• • # 

where X 0 , Y Q and Z Q are the velocity components, along X Q , Y Q 
and Z Q axes respectively. Expanding eqn. (3.27) gives the 
following three first order nonlinear differential equations : 

• 

X Q * V x cos© cos* + Vy(0cos* sin© - sin* )+V z (cos* sin©+0sin* ) 

(3.28) 

Y q * cos© sin* + Vy(0sin* sin© + cos* )+V z (sin* sin©-0cos* ) 

(3.29) 

Z - -V y sin© + V v 0 cos© + V cos© (3.30) 

O A. I Z 

Equations (3.28) - (3.30) after linearization and non- 
dimensionalization reduce to 

Dx q = U cos© cos* - up sin* + Ua cos* sin© , (3.31) 

Dy 0 = U cos© sin* + Up cos* + Ua sin* sin© , (3.32) 

and 

Dz q - -U sin© + Ua cos© . (3.33) 

X Y q Z Q 

where x Q = |-~ , y 0 * and z o ~ U~ * 



3.7 SUMMARY 


The trajectory is completely described with the help of 
nine differential equations (3.21) - (3.26) and (3.31) - (3.33) 
The six of these differential equations are of first order and 
the other three of second order. 12 initial conditions are 
required to solve these equations. 

They are at t + * o 
U ** 1 , 

a = p = 'Jf = 9 = 0= 'i = 0 = 0 = 0, 

and * y^ = z„ = 0 • 

o 1 o o 

The variation in atmospheric properties and wind which affect 
the trajectory significantly have not been considered so far. 
The corrections due to prevailing meteorological conditions are 
discussed in the next chapter. 



CHAPTER 4 


ATMOSPHERE AND ITS . EFFECTS 

4*1 ATMOSPHERIC STRUCTURE AND EFFECT OF COMPRESSIBILITY 

The aerodynamic derivatives appearing in eqns. (3.21) - 
(3.26) which depend on the Mach number would change in the 
course of flight of the projectile due to change in speed and 
also due to variation in temperature of the atmosphere with 
height. Further, the parameters involving density also change 
with height. In order to obtain the actual trajectory these 
effects should be suitably incorporated by considering proper 
variations in density and temperature with altitude and 
accounting for the compressibility effect. Specifications of 
ICAO atmosphere are assumed in the present study and the 
corresponding variations in density, temperature and speed of 
sound for an altitude upto 11000 m are given by 

p = p Q ( 1-0 .2255769 5xlO~ 4 h) 4 * 25587 , (4.1) 

T = T Q - 0.0065 H , ( 4 - 2 ) 

a = V YRT (4.3) 

where p is the density of air 

p is the density of air at mean sea level = 1.225 kg/m 
o 

T is the temperature of air 
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T c is the temperature of air at mean sea level 

= 288.15°K 

H is the altitude in metres 
a is the speed of sound in m/sec 

y is thb ratio of specific heats for air = 1.4 

R is the gas constant for air = 287.26 m 2 sec~^ °K~ 3 ‘ 


A similarity rule for axisymmetric bodies in the subsonic 
regime is used to account for the Mach number effect on the 
aerodynamic derivatives. According to similarity rule the 
pressure coefficient at any given point on the body, C ^ at a 
Mach number M a ^,is related to the pressure coefficient C \ ^ at 
another Mach number M a2 -by[2] 


„ 1 - ^2 
Pl = r^§i C ? 2 


(4.4) 


The aerodynamic derivatives that arise from the normal 
pressure distribution can be corrected using this factor. 
However, the drag coefficient which is mainly due to the tan- 
gential pressure distribution can not be treated similarly. 

In a finned-projectile both the fins and the body contribute 
to drag coefficient. The total drag coefficient Cq is given by 


C D * + ^ C D^W ’ 


(4.5) 


where 

(C D ) B is drag coefficient due to body, and (C Q ) W is drag 
coefficient due to fins . 
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The body contribution to the drag coefficient is given by 


( C D^B ~ C DP + C DF + C DB * ( 4 *6) 

where 

C DP 3 P ressu ^ e drag coefficient, 

Cgp = skin friction drag coefficient, 

and , 

C DB = base drag coefficient. 


The pressure drag coefficient C^p is small compared to 
skin friction drag and base drag for the type of bodies under 
consideration and hence is neglected. The other two components 
are obtained from the following relationships [3] : 


The skin friction drag coefficient Cpp is given by 

r - r 4.96 x 1(T 4 SL 0.032 , 1 0,1 S wet 

DF ( 1+0 .03Mg ) ^wet ^/i+o.i2Mg Re0 * 145 Swet Smid 

(4.7) 

where 


» porition of the wetted area over which laminar flow 
exists, 

S . = total wetted area, 
wet 

V L 

Re = Reynolds number = ~ , where v is the kinematic 

viscosity of air and V and L are as already defined. 


S mid 


« largest cross-section area, and 


M*- flight Mach number. 
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In eqn. (4.7) the contributions to skin friction due to both 
the laminar and the turbulent flow have been included. The 
first part is due to laminar flow and the second part is due to 
turbulent flow. The base drag coefficient is expressed in 
terms of C^p as 


r - 0*029 r B 
DB fC^ {d mid 


) 3 , 


(4.8) 


where 


d Q = base diameter, and 


d mid 36 ^ ar 9 es ^ diameter. 

The fin contribution to drag coefficient is given by [4] 

S 


(Cp) w = C f [l+ L( t/c)+100( t/c) 4 ] R^ s 


wet 
S ref ’ 


(4.9) 


where 


= 


T/c 

Z 


wet 

S ref 


turbulent flat-plate skin friction coefficient as a 
function of Mach number and Reynolds number, 

average stream wise thickness ratio of the fin, 

airfoil thickness location parameter, L = 2.0 for 

(t/c) max located at x^. < 0.30 c (x^ is the chord- 

wise position of maximum thickness) (4. 1.5.1) [4], 

wetted area of the fin, 

reference area, and 

lifting surface correction factor. 
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The lifting surface correction factor R^g is a function 
of Mach number. For this type of projectiles an approximate 
relationship between R^g and Mach number is given by 

R ls = 1,5 Ma- (4.10) 

The aerodynamic coefficients except the drag coefficient 
have been estimated for a Mach number of 0.7 using DATCOM 
techniques and the details are presented in Appendix B. They 
have been subsequently corrected for actual Mach number using 
the similarity rule at every point on the trajectory during the 
computation. The drag coefficient on the other hand has been 
estimated at every point on the trajectory using the method 
discussed above for the actual Mach number. 

4.3 EFFECT OF WIND 

Wind has a considerable effect on the trajectory of the 
projectile. In order to analyse this effect on axis system 
X Q ’ Y q ' Z q * parallel to X 0 Y 0 Z Q and moving with wind speed 
in the direction of wind is considered. The dynamical equa- 
tions (3.21) - (3.26) and the kinematic equations (3.31) - 
(3.33) derived in Chapter 3 are now valid in ' Y q ' axis 

system. To obtain the trajectory in X 0 Y Q Z 0 the effect of 

translation of X ’ Y * ,Z ' is accounted by including the wind 

o o o 

velocity components along X Q# Y Q and Z Q axes. If W, the wind 

vector has components W Y > Wy , W 7 along X , Y , Z Q axes then 

A o o o 

the kinematic equations (3.31) - (3.33) accordingly get 
modified to 
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Dx q = U cos© cos* - Up sin * + Ua cos* sine + r (4.11) 

v o 

Wy 

Dy 0 = U cos© sin* + U(3 cos* + Ua sin* sin© + , (4.12) 


z 

Dz q = -U sin© + Ua cos© + . (4.13) 

J v o 

If U? has components W , Vf , W along X,Y,Z axes then the angle 

A Y A 

of attack, angle of side slip and the relative air speed Vj^ are 
given by the following modified expressions : 


sin a 


sin p 



(4.14) 

(4.15) 



(4.16) 


Since angles a and p are small eqns. (4.14) and (4.15) can be 
approximated by 


and 



(4.17) 

(4.18) 


It can be seen from eqns. (4.17) and (4.18) that even 
though Vy = V z = 0 at t + = 0, angles a and p are not zero. 

This calls for a modification in the initial conditions which 
arise from wind. The initial values of these two angles are 
therefore, 
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and 






(4.19) 


(4.20) 


(4.21) 


In Fig. 4.1 the relative orientations of . X 0 Y C Z 0 axis system 
and W in horizontal plane are shown. The wind direction and 
line of fire are always referred to North which is also shown 
in this figure. If © w and © L are the angles between North and 
wind direction and North and line of fire respectively then 

W Xo> "yo’ W Zo are 9iven by 

W x = cos(e w - 0[_) COS0 O , (4.22) 

o R 

w Y(> “ UjJ sin (® W ' e L 5 ’ <4 ' 23) 

W_ = rr- cos(© - © L ) sin© . (4.24) 

z. 0 v R 


and 
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CHAPTER 5 

SOLUTION OF EQUATIONS 

5.1 GENERAL 

The eqns. (3.21) ~ (3.26) and (4.11) - (4.13) have been 
solved numerically on Dec-10 computer using the library 
routines written by the ’Numerical Algorithms Group' (NAG) [5]. 
Three methods, namely the Runge Kutta Mcrson method, the Adams 
method and the Gear method were tried and the Gear method was 
chosen based on time consideration. This routine integrates 
a system of first order, ordinary nonlinear differential equa- 
tions with suitable initial conditions using a variable -order, 
variable-step Gear method and returns the solutions at points 
specified by the user [5]. The accuracy of integration and 
interpolation are controlled by the parameters TOL and IRELAB 
which are discussed in the subsequent paragraph. The stiffness 
of equations, if any, is also accounted by this routine. 

4.2 ERROR CONTROLLING PARAMETERS 

The parameter TOL specifies the accuracy tolerance. 

The step size is automatically adjusted by the routine to 

-5 

meet the tolerance requirements. Values of TOL between lO 

—7 

and 10 were tried and it was noticed that the results were 

—5 

uneffected upto the 5th significant digit. A value of 10 
was chosen to keep the execution time small. 
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The other parameter which controls the error is IRELAB. 

It controls the error based on the value specified for TOL. 

At each step in the numerical solution an estimate of the 
local error, EST is made, For the current step to be 
accepted the following conditions must be satisfied : 

i N 9 

IRELAB =0 EST < TOL x max[l.O, w 2 Y(l) ] ; 

1M I=1 

IRELAB =1 EST < TOL % 

i N 9 

IRELAB =2 EST < TOL x max [esp, ^ E Y(l) z ] . 

N 1=1 

where eps is a small machine dependent number, Y is the matrix 
that contains all the variables of the set of differential 
equations being solved. 

The program requires on an average about 5-6 seconds of 
CPU time to compute a complete trajectory and write all the 
variables at steps of 500 non dimensional time. The time 
required for writing the results at all intermediate points 
which is about 30 - 40 °/, of the total time can be eliminated if 
only the results at the end point are required. The program 
in its present form requires 70 blocks on the disc. 

4.3 PHYSICAL PARAMETERS OF THE TRAJECTORY 

Besides the variables appearing in the equations of 
motion, the parameters of interest along the trajectory are 
the range, height and the angle of graze at the end point. 
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all referred to the horizontal plane passing through the 
point of launch. 

Expressions for range and height are obtained referring 
to Fig. 4.1 and are given by 

Range = L(x 2 cos 2 © 0 +y 2 +z 2 sin 2 e o +x Q z o sin 2© 0 ) 1//2 , (5.1) 

Height^ L(x rt sin 9 - z^ cos© A ).. (5.2) 

^ o opo 

The angle of graze a is given by 

6 = 2 “ cos (-cos© cos $ sin© o -sin0 cos© 0 ) (5.3) 

In the present study the equations have been solved for a 
120 mm artillery projectile and the results are presented in 
the next chapter. 
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CHAPTER 6 

RESULTS AND DISCUSSION 


6.1 GENERAL 

Time histories of angle of attack, angle of side slip, 
angle of roll, angle of pitch, angle of yaw and the lateral 
deviation and also the variations of range, height and lateral 
deviation for different angles of projection, muzzle velocities 
and wind conditions have been studied. The angle of graze at 
the terminal point of the projectile has also been determined 
for different angles of projection. Due to uncertainties in 
the estimated value of the aerodynamic coefficients, a para- 
metric study has been carried out to ascertain the sensitivity 
of the trajectory to variations in these coefficients. The 
results are presented in Figs. 6.1 - 6.29 and discussed in the 
subsequent paragraphs. 

6.2 RANGE AND HEIGHT 

Figures 6.1 and 6.2 show the trajectories for several 
muzzle velocities and angles of projection. The range table 
values of the range and height are plotted for comparison. For 
60° projection angle the trajectory of an equivalent point 
mass is also presented. Percentage deviations in the overall 
trajectory parameters compared to the range table data are 
presented in Table 6.1. 
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It is seen clearly from Fig. 6.1 and Table 6.1 that the 
finite size of the projectile and its associated aerodynamics 
has a significant contribution for the trajectory. It can also 
be seen that the computed results are very close to the range 
table results. Range and height increase with increase in 
muzzle velocity where as with increase of angle of projection 
the range decreases and height increases. 

In Fig. 6.3 variations of range and height versus muzzle 
velocity are presented. They are also compared with the range 
table data. The plots show very close agreement of both range 
and height with the present computed results for all muzzle 
velocities. It is further noted that while the estimated range 
is lower than the range table data for all muzzle velocities, 
the estimated height is lower than the range table for muzzle 
velocities below 205 m/sec and higher than the range table data 
for muzzle velocities above this. 

In Fig. 6.4 variations of range and height with angle of 
projection are shown. For comparison the range table data 
are also plotted and it can be observed that the computed 
results for both range and height are very close to the range 
table results. 

6.3 SPEED, ANGLE OF PITCH AND ANGLE OF ATTACK 

In Figs. 6.5, 6.6 and 6.7 the time histories of non- 
dimensional speed, angle of pitch and angle of attack are 



Table 6.1 Comparison of computed results with range 
Table Data 
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presented. The range table value of the terminal velocity at 
the point of graze is also shown. It is observed that the 
terminal velocity from computation agrees very closely with 
the range table value. The minimum speed of the projectile 
along the trajectory occurs at the vertex, i.e., the point 
where height is maximum. From Fig. 6.6 it can be seen that 
the angle of pitch remains negative through out the trajectory, 
of the projectile. It builds up slowly in the initial part 
up to about ^-th of the trajectory and thereafter its rate of 
build up increases till towards the end of the trajectory 
where once again the build up rate becomes small. Angle of 
attack variation in Fig. 6.7 shows that it' is oscillatory and 
the oscillations get damped out in about ^th of the trajectory, 
after which it practically remains at zero. Since 0 variation 
presented in Fig. 6.6 does not indicate this type of oscillatory 
nature, these transient oscillations in a are due to the heav- 
ing motion of the projectile. 

6.4 EFFECT OF WIND 

6.4.1 Effect of Range Wind 

6.4. 1.1 Range 

Effect of range wind on range is shown in Fig. 6.8 and 
the results are compared with range table data. It is noticed 
that while a tail wind increases the range a head wind results 
in a decrease of it. In either case the variation in range is 
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linear with wind speed. Comparison with the range table 
results shows that the present method slightly underestimates 
the range for head wind and the difference reduces as the wind 
speed increases. As for the effect of tail wind, the 
present study gives an underestimate of the range for wind 
speeds upto about 10 m/sec and for speeds above this it gives 
an overestimate of the range. 

6.4. 1.2 Angle of Attack 

Figure 6.9 brings out the effect of range wind on angle 
of attack. It is seen that the angle of attack variation is 
characterised by high frequency transient oscillations. The 
peak amplitude of these oscillations is highest for head wind 
and lowest for zero wind as can be seen from Figs. 6.7 and 
6.9. It is also noticed from these figures that these tran- 
sient oscillations subside earlier for tail wind than for 
the head wind. 

6.4.2 Effect of Cross Wind 

Throughout the study a right cross wind is considered. 

For a left cross wind the results will be mirror images of 
those for right cross wind. 

6. 4. 2.1 Lateral Deviation 

In Figs. 6.10 and 6.11 the lateral deviation Y Q is 
plotted as a function of non dimensional time and the 
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variation of the terminal value of Y Q with cross wind velocity 
is presented in Fig. 6.12* It can be seen from Figs. 6.10 and 
6.11 that build up in is quite slow and gradual initially. 

For lower wind velocities the trend continues to be gradual 
whereas for higher wind velocities the build up in the later 
part of the trajectory becomes quite rapid. It even seems 
to reach a maximum value after which it starts decreasing as 
can be noticed from the curve for V = 248.1 m/sec and 9 rt = 70° 
in Fig. 6.11. In order to understand this type of behaviour 
of lateral deviation, the two major contribution for Y Q have 
also been plotted for 9 Q - 70° on the same figure. The wind 
contribution increases linearly arid is positive for right cross 
wind throughout the trajectory whereas the contribution due to 
yaw is negative and increases linearly only in the first half 
of the trajectory. In the second half of the trajectory this 
term decreases rapidly and starts increasing again after some 
time. Because of this behaviour of the contribution due to yaw, 
Y Q increases slowly and almost linearly in the first half and 
rapidly in the later part. This suggests that the variation of 
terminal Y Q with wind velocity may not be monotonic for higher 
wind velocities and large angles of projection and Fig. 6.12 
indeed, confirms this. In this figure the range table results 
are also plotted. For angles of projection upto 60°, Y Q varia- 
tion is almost linear with wind speed upto about 10 m/sec 



whereas for © 0 > 60° this monotonic variation breaks down. 

This is because of the yaw angle component as discussed 
earlier. In range table the wind effect has been compiled 
based on the concept of equivalent constant wind. This 
implies a linear dependence of Y Q on wind velocity and does 
not take into consideration the contribution of yaw to lateral 
deviation. Since this second contribution is quite significant 
at large angles of projection and wind speeds as has been 
noted, the lateral deviation as compiled in range table for 
these conditions can be highly unreliable. 

6. 4. 2. 2 Angle of Sideslip 

Figure 6.13 shows the time history of angle of sideslip p 
with cross wind. The angle of sideslip is nonzero initially 
due to cross wind and it damps out quite rapidly to zero steady 
state value after going through transient oscillations as seen 
from the figure. This indicates that the projectile aligns 
itself to the resultant wind direction quickly after going 
through lateral oscillations initially. 

6. 4. 2. 3 Angle of Roll and Angle of Yaw 

Time histories of angle of roll and angle of yaw are 
presented in Figs. 6.14 and 6,15 for different angles of 
projection and muzzle velocities. From Fig. 6.14 it is seen 
that for 9 = 45° and V » 153.0, m/sec, both 0 and ¥ show 

oscillatory behaviour initially. The oscillation die out in 
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about ^th of the trajectory after which they start slowly 
increasing. The magnitude of 0 is small compared to that of S 
(0 is very close to zero) initially. Towards the end of the 
trajectory both of them build up rapidly and reach maximum 
values. Beyond this they start reducing equally rapidly 
and seem to cross the time axis if the trajectory had continued. 
This change of sign of 0 and ¥ is realised for higher muzzle 
velocities and angles of projection as can be seen from Fig. 
6.15. In this figure also t initially the magnitude of ¥ is 
much larger than that of 0 . This fast build up in the magnitude 
of 4 initially is because of yawing of the projectile in res- 
ponse to nonzero side slip by virtue of its directional stabi- 
lity C . The projectile undergoes rapid yaw oscillations till 

P 

it aligns to zero sideslip. The corresponding initial oscilla- 
tions in p have been shown in Fig. 6.13 and have already been 
discussed. It is also noticed that over almost the entire 
trajectory 0 and ¥ maintain opposite signs. 

6. 4. 2. 4 Angular Velocity Components 

Figure 6.16 presents the time histories of the angular 
velocities p,q and r. The roll rate p shows a build up in 
oscillations initially which subsequently die out. Variation of 
q is negative and non oscillatory and reflects the variation 
of 0 in Fig. 6.6. The yaw rate r however, shows oscillations 
of much larger frequency than p and dies out much more rapidly 
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to zero than p. This is due to rapid oscillatory variations 
in ^ discussed earlier. In order to study the correspondence 
between the angular velocities and the time rate derivatives of 
the Eulerian angles their time histories are plotted together 
in Figs. 6.17 and 6.18. It can be seen from Fig* 6.17 that p 
and 0 variations are different except between 1500 and 6000 non- 
dimensional time. This can be understood clearly from the ex- 
pression for p in eqn. (2.2). In the initial part of the tra- 
jectory which is characterised with the presence of transients, 
the magnitudes of both the contributions for p namely, ¥ sin© 

and 0 are comparable and hence p and 0 vary differently. 

+* • 

After about t == 1500 while 9 is still small U 1 also becomes 

very small and the contribution to p due to ¥ sin© becomes 
relatively insignificant. Therefore, p and 0 show almost iden- 
tical variations in this part of the trajectory. Subsequently 
for t + greater than about 6000, 0 also builds up and results 
in different variations of p and 0 . The variations of r and ¥ 
in Fig. 6.18 are identical initially after which r becomes 
almost zero while $ continues with a small value and eventually 
becomes zero. 

6.5 ANGLE OF GRAZE 

In order to study the orientation of the projectile with 
respect to the horizontal plane at the point of graze, the 
variation of angle of graze 6 with angle of projection is 
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plotted in Fig. 6.19. The angle of graze variation with @ Q is 
linear throughout. If the angle of yaw is small at the terminal 
point, the magnitude of 6 is equal to |o| - © . Since 6 is 
found to vary linearly with 0 Q it can be concluded that (©| at 
the point of graze also varies linearly with © Q . The effect of 
wind on angle of graze is also shown in this figure. Increase 
in wind velocities above 4.0 m/sec did not show any noticeable 
change in the angle of graze. 

6.6 PARAMETRIC STUDY 

The variations in the trajectory parameters range, height, 
lateral deviation, angle of attack and angle of side slip due 
to variation in the aerodynamic coefficients are presented in 
Figs. 6.20 - 6.29. 

6.6.1 Effect of Drag Coefficient 

Figure 6.20 shows the sensitivity of range, height and 

lateral deviation to drag coefficient C v . It is noticed that 

A o 

both the range and height increase with decrease of drag 
coefficient and decrease with increase of it. The effect of 
drag coefficient on lateral deviation is seen to be insignifi- 
cant. Drag coefficient is the most significant parameter in 
determining the range. 

The effects of other aerodynamic coefficients namely, 

r c . C . C . C , C and C, on the range, height and 
* Z* 1 Z ’ m * m« r m l n 

a a q a a q P 
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lateral deviation are shown in Figs • 6*21 — 6*27 . It is seen 

that the variations of these coefficients have no significant 

effect on the range and hgight. However, the effect on lateral 

deviation is noticeable in some cases. Y decreases continuous— 

o . 

ly with increase of C z , C z and C m . As for the effect of C- , 

a q q a 

the lateral deviation is minimum at the base value of C-, and 

increases on either side of the base value as shown in Fig. 6.22. 


A similar variation of Y Q is seen in case of C, also (Fig. 6. 27). 

P 

With C , Y has its maximum at the base value and decreases for 
a 

other values of C on either side of it. The effect of C » is 

rn m* 

a a 

seen to be practically insigificant . However, at lower values 

of C the lateral deviation obtained is slightly more than 
m a 

that obtained at the base value. 


6.6.2 Effect of C on Angle of Attack 

m a 

Figure 6.28 shows the time histories of angle of attack 

for different values of C m . The amplitude and frequency of the 

m a 

transient short period oscillations increase with increase of 

C . The increase in frequency is due to increase in stiffness 
in, 

OL 

(in pitching motion) arising from C . It is also noticed that 

a 

for lower values of C the short period oscillations are 

m a 

followed by oscillations of relatively larger period which con- 
tinue over a large part of the trajectory and this behaviour 

is almost not there for higher G . 

m a 
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6.6.3 Effect of C on Angle of Side Slip 

"(3 

In Fig. 6.29 time histories of angle of side slip for 

different values of C are plotted. It can be seen that the 

0 

frequency of oscillations increases due to increase of C 

n P 

from -0.792 to -2.286. However, the curve for C = -1.585 

n p 

shows very few oscillations and is better damped. This is 

likely to be due to the fact that a variation in C changes 

P. 

the damping ratio and C = -1.585 may be closer to critical 

0 

damping than the other two. 
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CHAPTER 7 

CONCLUSIONS 

Obtaining the trajectory by solving equations of motion 
jis a generalized approach and the method can easily be applied 
and extended to obtain the trajectory of any projectile with 
similar or modified configurations. The results of the present 
computations have shown that the angle of attack a, angle of 
side slip (3, angle of roll 0 and the time derivatives $ , 0 and 
0 are very small and this justifies the assumptions made in 
simplifying the equations of motion. 

The results of the present study show a good agreement 
with the range table results for range and height. But the 
lateral deviations in the two cases differ considerably. The 
differences in the computed results and the range table data 
can be attributed to the following : 

. Uncertainties in the compilation of range table data 

which can be due to variations in the gun performance, 
limitations of measuring instruments and the errors in 
Interpolation. 

. The scheme used in the preparation of range tables to 
reduce the fired-data to the standard meteorological 
conditions itself is unsatisfactory and may introduce 
errors in the range table results* 
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. .The range table assumes a linear dependence of the 

lateral drift on cross wind and neglects the contribu- 
tion arising from the yaw angle term. The second 
contribution is quite significant and cannot be 
ignored. 

. Uncertainties in the estimation of the stability deriva- 
tives of the projectile in the present study. 

Use of similarity rule to account for the compressibility 
effects is not exact and hence might introduce errors in 
the aerodynamic derivatives. 

The parametric study carried out by varying the aerodynamic 
coefficients shows that the trajectory is significantly affected 
by only a few of these derivatives. It is found that the drag 
coefficient has the most dominant effect and the rest of the 
derivatives have negligible effect on the overall trajectory. 

It is therefore, imperative to determine the drag coefficient to 
a high degree of accuracy. 

In this study the ICAO standard atmosphere has been used 
but the method is applicable in any atmosphere whose structural 
properties are known. 

The speed and accuracy of computation of the present method 
strongly recommend its application in the modern environment of 
battle field. The program to solve the equations is small in 
and occupies only 70 blocks on the disc in the present 
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form. Since the NAG routines used to solve the differential 
equations are also available on small computers, it would be 
possible to implement this model on a microprocessor. A 
ruggedised model of such a microprocessor will further advance 
its suitability for application in the battle field. 
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CHAPTER 8 

SCOPE FOR FURTHER WORK 

The present method can be successfully applied to study 
many other aspects of the motion of a projectile which could 
not be carried out due to time restriction and non availability 
of sufficient aerodynamic data. It has a wide scope and a few 
suggestions in this regard are as follows : 

The present model is restricted to subsonic speeds. The 
technique can be extended to transonic and supersonic 
speeds by choosing proper aerodynamic coefficients 
obtained either analytically or from wind tunnel tests. 

. The equations of motion can be applied for rocket- 

assisted projectiles by including the thrust effect. 

. The effect of rotation and curvature of earth can be 

incorporated which are significant only for long range 
projectiles. 

. Probabilistic study on the point of fall can be carried 
out based on some given distribution of muzzle velocity, 
angle of projection etc. 

• The present model can be optimized, an error analysis 
made for smaller machines and finally implemented on a 


microprocessor 
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DERIVATION OF EQUATIONS OF MOTION 


% 


m 


The equations of motion presented in Chapter 3 are derived 

step by step here. Linearization and non dimensionalization of 

’ ■ 

these equations are also discussed in detail. 

\ 

X COMPONENT OF FORCE 


Equation (3.1) is resolved along X axis and transformed 


using eqn. (3.9) to get 

q ^V r ^Y = F AX + F GX * 


.. 


or 


ft (m V x } + qmV z -rmVy = + F, 


GX 


Substituting for V^, V y , V z from eqn. (2.1), for q,r from eqn. 
(2.2) and for F^ and Fq^ from eqns. (3.18) and (3.20) res- 
pectively It follows that 


m 


r r 2 2T dv V( sin20C if + Sin2p Ht } 1 

[ V 1-sin a-sin P ^t “ ?====?- J 


2/l-sin 2 a-sin 2 p 


+ m(©cos0 + i sin0 cos©) Vsina-m(* cos© cos0 - ©sin0)VsinP 


£ - mg(sin© 0 cos© cos* + cos© 0 sin©), 


or 


. 


§Jfii 
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m^fl-sin a-sin p dt ” 2 + m 9(sinQ o cos© cos! + co.s© Q sin©) 

mV( sin2a + sin2p 

= ‘ ‘T "2 2™* ~ m(©cos0 + ! sin0 cos©) Vsina 

Zyl-sm a-sm p 


• * 

+ m(! cos© cos0 - Qsin0) Vsinp 


(Al) 


Y COMPONENT OF FORCE 

Equation (3.1) is resolved along Y axis and transformed 
using eqn. (3.9). This gives 

M y + r ^)T P ^Z = F AY + F GY * 

or 

J 

dt ( mV y ) + rmV^ — pmV^ — ^AY ^GY * 

Substituting for ^x ,V Y* V Z from ec * n * for P» r from eqn. 

(2.2) and for F^y and F Q y from eqns. (3.18) and (3.20) res- 
pectively it follows that 

msinp ^ + mVcosP |£ + m( ! cos© cos0 - ©sin0) V /l~sin 2 a-sin 2 p 
-m($~ 4 sin©) Vsina = f V^Cy P+Cy §7 + C Y r f? ) 

+ mg[0(cos© Q cos©-sin© 0 sin© cos f )+sin© Q sin!], 
or 

msinp + mVcosp + mv/l— sin a— sin p ( 4 cos© cos0) 

- | ^(Cy P+Cy. + Cy^ f COS© cos0)-mg[0(cos© (EJ) cos© - 
sin© sin© cos!)+sin© 0 sin! ] = mVsina(0 -4 sin©) 

+ mvy l-sin 2 a-sin^p(©sin0) - 5 C y^ 2V ® sin ^ 
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Z COMPONENT OF FORCE 


Equation (3.1) is resolved along Z axis and transformed 
using eqn. (3.9) to obtain 

Mz + P^y*“ qM - p + p 

x. Mi r AZ -t* r GZ , 


or 


dt ^ mV Z^ + praV Y ~ qmV x = F az + F qz . 

Substituting for V x ,Vy,V x from eqn. (2.1), for p,q from eqn. 

(2.2) and for and Fqz from eqns. (3.18) and (3.20) res- 
pectively results in 

msina ^ + mVcosa ^ + mVs inp((ZS - « sin©) 


•m(©cos0 + ¥ sin0 cos©) V^l-sin 2 a-sin 2 /3 = | V 2 ?^' a + 

a 


T 

C 7- 5v + C Z 2v (® COs $ + ? sin0 cos©) ]+mg[cos© cos© - 
a q 

sin© (cos 'Si sin© + 0sin ¥ )], 


or 


msina + mVcosa mv/l-sin 2 a-sin 2 p(©cos0) 


Ht 

I V 2 S( C z a+c z . ^ ^ ° 


-a 


+ C Z 5v ecos0)- m g[c os ©o COS© 

q 


sin© ( cos f sinQ+0sin f ) ] = -mVsinp(0 - 5 sin©) + 


m 


\pj 1 — sin 2 a-s in 2 p ( ^ sin 0 cos©) + 5 V 2 ^ C z 7 ^ ^ sin0 cos© 

q (A3) 
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X COMPONENT OF MOMENT 

Equation (3.2) is resolved along X axis and transformed 
using eqn. (3.9) to get 

H x + qH^ - rHy = M x . 

Substituting for H x , Hy and H z from eqn. (3.7) gives 

j 

dFt "^* * M X * 

Substituting = ly and for p and M x from eqns. (2.2) and 
(3.18) respectively it follows that 

I x (3 - J sine) * i V^SL C x ~ (i - i sine), 

or 

X x - I x sine *j-| - | ^SL 27 (H - sinS I?) 

= i x cosef||. U4, 

Y COMPONENT OF MOMENT 

Equation (3.2) is resolved along Y axis and transformed 
using eqn. (3.9). This gives 

By + rH x - pH^ — My • 

Substituting for H x ,Hy,H z from eqn. (3.7), for p,q,r from 
eqn. (2.2) and for My from eqn. (3.18) results in 
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I Y dt (9cos ^ ^ c °s6 sin0) + (S-fsin0)(i cos© cos0-©*sin0) d x “I z ) 

= I V 2 SL [C a +C m . ~~ + c m |y (9COS0+1 cos© sin0)], 
oc a q 

or 

! y cos0 - | V^SL(C a + C § + C ^ Ocos0) 
u v. a a q 

= IySin0 ~ IySin0 cos© - Iy ~j— • (cos© cos0 

- sin0 sin© -j^) + — V^SLC^ Ijy ¥ sin0 cos© - 

q . 

(0-Wsin©)^ cos© cos0 - ©sin0) d x ~I z ) (A5) 

2 COMPONENT OF MOMENT 

Equation (3..2) is resolved along Z axis and transformed 
using eqn. (3.9) to get 


+ pHy - qH^ = M z 

Substituting for H x ,Hy,H z from eqn. (3.7), for p,q,r from eqn. 
(2.2) and for M z from eqn. (3.18) it follows that 

I 4r-( i cos© cos0-©sin0) + ($- ^ sin©) (©cos0+ S cos© sin0)(Iy-I x ) 

Z. at 

= £ V^LtC p+C . §n + C (i cosS cos0-Qsin0) ] , 

^ *p **p ^ x 
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or 

2 ' # 

I z cos© COS0 0 - f ^SUC^P+C^ §£f + ^ i cose cos 0) 

= -(0- ^ sin©) (9cos0+ ¥ cos© sin0) (I y -Iv)+I z sin0 

dt z 

+ cos0 "5? + *z HT ( C0s9 sin ^ + cos ^ sin ® H?) 

- | V^L C ©sin0 ( A6) 

r 

LINEARIZATION 

Equations (Al) - (A6) represent six scalar equations of 
motion of the projectile and are arranged in such a way that 
most of the higher order terms appear on the right hand side ! 

of the equations. All the first order terms have been written 
on the left hand side. Since the magnitude of the second and 
higher order aerodynamic terms is insignificant relative to the 

i 

first order aerodynamic terms, all the second and higher order 

I 

terms on the right hand side of eqns. (Al) - (A6) are neglected, j 
The equations are further linearized by making the following 
simplifications : j 


sin0 =0 ; 

cos0 = 1 - 2"“ 

sina = a ; 

„2 

t ft 

cosa - 1 - 2 ~ 


B 2 

sinp * p » 

cos (3 =1-5- 
= i - (2^-) .- 

Vl-sin 2 a-sin 2 p 
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After making these substitutions the second order terms in 
angles 0, cc and p are also neglected. Thus eqns. (Al) - (A6) 
can be simplified to get 

m Ht ~ "2 S( + mg(sin© Q cosQ cosf + cos© Q sin©) = O, (A7) 
mp + mV || + mV cos© H 5v + 


Cy 5v cos ® -mg[0(cos© Q cos© - sin© 0 sin© cos's ) + 

Y r 


sin© 0 sin* ] = 0 


(A8) 


L 

na + mV |f - mV || - § As( C^a+ 0 Z . If + c Zq - 2 y at 


) 


- mg[cos© 0 cos© - sin0 Q ( cos « sin© + tfsin'S )] = 0 , (A9) 

l y it - l x sin 0 A- - § As LC X & (H - sino ff) = 0 . ( AlO) 


dt 

* . 

t d 2 © £ v^St (C a + C ^7 = 0 * 

X Y ^2 - 5 ^ SL ^m/ 2V dt m q 2V dt 


(All) 


_ d 2 * P V 2 QT / r o + c -i.^4+C 4b cos © H+) = 0 
1 ,, 7 cos© — - -wr - k V SL vC P + L n * 2V dt n 2V at 

Z dt 2 “ P P (A12) 



To non-dimens ionalize eqns. (A7) 
time factor is defined as 


- (A12), a projectile 


(A13) 
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where, 

t : projectile time factor 
L : characteristic length 
V Q : initial velocity of projection . 
The non dimensional time t + is given by 


or 



Differentiating eqn. (A14) with respect to t gives 



(A14) 


(A16) 


The time derivatives in eqns. (A7) - (A12) are now replaced 
by the derivatives with respect to non dimensional time t . 
Using concise notations for the coefficients of the variables, 
the eqns. (A7) - (A12) with non dimensional time as independent 

variable are finally expressed as follows : 


DU - C x -U 2 + F( sin@ 0 cosQ cos* + cos© 0 sin©) - 0, (A17) 

2 . 

u- k + k ( VV U P 

-F(sin© 0 cos© cos* + cos© 0 sin©)P - F[0(cos© o cos© 
sin© Q sin© cos* ) + sin© 0 sm * ] - 0 » 


(A18) 
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u - k C Z-) UD “-( 1+ k C 2 >UD ® + k (C X - C Z a )u2a 

- F(sin© 0 cos© cosf + cos@ o sin©)a - F[cos© 0 cos© - 


s in© Q ( cos ¥ sin© + 0sin ¥ ) ] =0 


D 2 0-sin© D Z % -K-, C, UD0 + K x C, sin© UDi = 0 

P P 

D 2 e-S C m . UDO-I^ C m UD9 - 2K 2 C U 2 a = 0 
a q a 


(A19) 

(A20) 

(A21) 


cose D 2 * C n . UDP - K2 C n cosS UD* - 21C, C U f> = 0 


(A22) 


where, 


D = ~ * 
dt + 


U = 


V 

^ ' 


m 


^ = PSL ' 


3^ t 

V o 


Ki = 


PSL _ 
“4T x ' 


PSL 3 PSL 3 

K 2 =^|fc= -ir 
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APPENDIX B 

PROJECTILE BODY PARAMETERS AND ESTIMATED AERODYNAMIC 

COEFFICIENTS 


BODY PARAMETERS 


The results have been computed in the present study for a 
120 rum projectile which has the following body parameters ; 


Total length : 0.6746 m 

Location of centre of gravity from base ; 0.4037 m 

Maximum diameter : 0.1196 m 

Base diameter : 0.0455 m 

Length of the fin : 0.0615 m 

Width of the fin : 0.042 m 

Number' of fins : 12 

Weight of the : 13.216 Kg 

projectile 


Moments of Inertia 



I 


Y 



0.02580 Kg m 2 
0.24610 Kg m 2 


AERODYNAMIC COEFFICIE NTS 

The aerodynamic coefficients have been estimated for the 


following conditions : 
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Characteristic length 
Characteristic area 
Speed 


5 0.675 m 

: 0.011122 m 2 
: 236.4 m/sec 


The estimated values of the aerodynamic coefficients are as 


follows : 

C = -0.151 [3], (4. 1.5.1) [4] 

C_ = -8.44 (4.3. 1.2) [4] 

z a 

C„ = -0.723 (7. 3. 4.1) [4] 

Z a 

C„ = -9.12 (7. 3. 1.1) [4] 

Z 

q 

C, = -0.6 (7. 3. 2.2) [4] 

l p 

C = -1.585 (4. 2. 2.1) [4] 

% 

C = 0.047 (7. 2. 2. 2) [4] 

m* . . 

a 

c s -2.754 (7. 3. 1.2) [4] 
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FIG* 8*10 LATERAL DEVIATION VS NON DIMENSIONAL TIME) %~ 1£*8 1 M/SEC, 
W - 4*0 M/SEC # CROSS WIND* 



NOliVlAia 






iliiss 

HI 







FIG-6. 15 VARIATION OF ANGLES OF ROLL AND YAW WITH NON DIMENSIONAL TIME ; 
So *60 DEG, Uo* 248-1 M/SEC, W*4-0 M/SEC, CROSS WIND- 








FIG. 6-10 VARIATION OF T AND V WITH NON DIMENSIONAL TIM£;0o-6ODEG 
VU *248-1 M/SEC» Wm4-0M/SEC, CROSS WIND- 
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FIG.6-20 EFFECT OF DRAG COEFFICIENT ON RANGE, HEIGHT ANO LATERAL 
DEVIATION •, 0o=6O DEG , 248 - 1 M/SEC, W» A'O M/SEC , CROSS WINO' 
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FIG- 6*28 EFFECT OF ON TINE HISTORY OF ANGLE OF ATTACK; 0f8OOE6 

-248*1 M/SEC, 






FIG. 6- 2 9 EFFECT OF Cn ? ON TIME HISTORY OF ANGLE OF SIDESLIP; 
*248-1 M /SEC, W. 4-0 M/SEC, CROSS WIND. 



